New Model and Numerical Test of a High Temperature Pairing Mechanism in Stripes 
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We introduce a new model and mechanism of high temperature pairing in stripes. We propose 
a way to unambiguously test it by numerical simulations. For example, the implementation of our 
mechanism in a 6-leg t-J ladder model has the effect of making the spin gap of the doped 6-leg 
ladder be about 4 times bigger than the spin gap of the same ladder at half filling. 



The cuperate superconductors manifest a remarkable 
variety of unusual and unexplained phenomena. Let us 
focus on superconductivity. As highlighted in their name, 
the high temperature superconductors are distinguished 
by a high pairing energy scale ~ SbmeV. The most 
significant missing element in our understanding of super- 
conductivity in the cuperates is the microscopic pairing 
mechanism. Known microscopic characteristics, such as 
the strength of spin exchange interactions, J w 1500K, 
restrict candidate mechanisms. 

Therefore, from a theoretical point of view we ask what 
models can, in principle, produce such high pairing en- 
ergies over the range of doping as observed in the cuper- 
ates? Our aim in this paper is to provide such a model 
based on stripes. We also propose a straight foreword 
numerical test of our proposed mechanism. The model, 
and the main effect which is in operation, is distinct from 
previous proposals. 

The paper is organized as follows: First, we present our 
model and propose a numerical test to verify or negate 
the operation of the core effect at work (we call it an "en- 
vironment spin decoupling effect" ) . Then, we use a two- 
chain model to illustrate the environment spin decoupling 
effect using tractable analytical methods. We conclude 
with a discussion of the implications of our model within 
the general phenomenology of stripes and superconduc- 
tivity in the cuperates. 



STRIPE MODEL AND NUMERICAL TEST 

Experimentally, the spin gap in cuperate superconduc- 
tors is nearly independent of doping 6 in the range 0.06 < 
(5 < 0.15 (with a magnitude = 2 A (0) 800K ~ J/2 
Q). In the same doping range, the distance between 
stripes changes by a factor two. Therefore, any mech- 
anism of pairing via stripes must be independent of the 
distance between stripes. The difficulty of achieving such 
a mechanism can be appreciated by thinking of stripes 
as a periodic arrangement of ladders: 

The mechanism of pairing in stripes is supposed to op- 
erate on each stripe independently. Therefore, the same 
mechanism should manifest itself in a single ladder model 
where the number of legs of the ladder is determined by 
the stripe period. Since the stripe period change from 
4 to 8 lattice spacings between doping 6 — 0.12 and 



FIG. 1: Schematic picture of a stripe state. Darker shaded 
regions represent higher hole consentration. Arrows repre- 
sent magnitude of spin moments on Cu sites. Periodic charge 
stripes (hole rich) are separated by nearly undoped regions. 



6 = 0.06, it implies that a mechanism of pairing in stripes 
should be such that the pairing of mobile fermions in the 
corresponding 4- leg and 8-leg ladder models be the same!! 
Indeed, the model and mechanism which we propose is 
aiming to achieve exactly that. 

In order to have a sharp goal in mind, we pose the fol- 
lowing theoretical question: Is it possible to have a spin 
gap in an even-leg doped ladder which is several times 
bigger than the spin gap of the same undoped ladder? 
Numerical investigation '3| seem to indicate that the an- 
swer is negative on uniform ladders. In this paper, we 
investigate non- uniform ladder models. 

One idea in this direction is to make the bond interac- 
tions inhomogeneous. Arrigoni & Kivelson considered 
a 4-leg Hubbard ladder model where the bare hopping 
interaction was dimerized transverse to the ladder (while 
maintaining the same uniform hopping interactions along 
all the legs). They showed that indeed the maximum spin 
gap of the doped ladder may increase substantially. 

We introduce an inhomogeniety not in the interactions 
but instead in the on-site chemical potential, e.g., for 
t-J type ladders The model is characterized by three pa- 
rameters {t, J, /i}. In Fig. 2 we draw a 3- leg and 6-leg 
ladder versions of our models. To begin with, we shall 
discuss the ladder model from a purely technical perspec- 
tive and defer the discussion of its association with the 
stripe state to the concluding section of this paper. The 
simplest anisotropic ladder model consists of a chemical 
potential shift /i in two neighboring legs of an N-leg lad- 
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FIG. 2: The empty circles denote some incommensurate elec- 
tron filling fraction. The other legs are half-filled or nearly 
half-filled and manifest significant local antiferromagnetic cor- 
relations. In this figure we draw the examples of anisotropic 
3-leg and 6-leg ladders, and the associated phase diagram at 
moderate doping. 



der. Consequently, holes are preferentially concentrated 
in these two legs. Hence, upon doping by holes away from 
half filling, the chemical potential difference leads to an 
inhomogeneous distribution of the holes. We argue that a 
significant renormalization of interactions leads to an ef- 
fective phase separation of the system into two incoherent 
subsystems at intermediate temperatures T2 < T < T* 
- One is a doped 2-leg t-J ladder, and the second is an 
undoped {N — 2) leg ladder. This decoupling effect is the 
main effect which leads to a high energy pairing scale. It 
is a direct consequence of the inhomogeniety of the hole 
distribution in the ladder due to the confining potential 
whose strength is characterized by fj,. At lower temper- 
atures T < T2 coherence is established by the proximity 
effect. As we discuss in section- IV, this lower tempera- 
ture coherence is related to the superconducting Tc- 

First we focus on zero temperature properties which 
can be checked by numerical simulations and serve as a 
qualitative and quantitative test of our claimed princi- 
ple mechanism. Consider the example of a 6-leg ladder 
(see Fig. 2). Upon doping, the system separates into a 
doped 2-leg ladder and an undoped 4- leg ladder, i.e., 
the minimal spin gap in the system is increased to the 
order of a 4- leg ladder spin gap w J/6 (i.e., a five 
fold increase compared with the undoped 6-leg ladder gap 
As ~ J/20!). This can be check by numerical simulation 



methods such as DMRG and CORE An additional 
conjectured property of our model is the following: The 
pairing energy of the holes is that of a two leg ladder, 
and largely independent of the spin correlations in the 
remaining undoped legs of the ladder. This should also 
be amenable to numerical experiments. 



ID MODEL WITH COMPETING INTRA- BAND 
AND INTER- BAND INTERACTIONS 

The above proposed numerical experiments can quan- 
titatively test our proposed model, yet, in the following 
section we would like to develop the reader's analytical 
intuition about the competition and renormalization of 
interactions in the model. Instead of directly thinking 
of the dynamic decoupling of a ladder system into two 
incoherent sub-systems, let us think of the reverse route: 
Consider two ID systems - one is a doped 2-leg Hubbard 
ladder and the second is an undoped or lightly doped N- 
leg Hubbard ladder. Turning on hopping interaction t± 
between the two ladders, there is a competition between 
the intra-ladder interaction which maintains the spin gap 
(and pairing) on the 2-leg ladder and inter-ladder inter- 
actions which form coherence across the whole coupled 
system. Similar issues of confinement/deconfinment have 
risen in the context of coupling equivalent ladder systems 
where the tuning parameter is inter-ladder interactions 
(e.g., t± ^ t 01 J± ^ J). In contrast, our model starts 
with uniform bare interactions (tj_ ^ t, J± ~ J) while 
inhomogeniety of doping induced by tuning an inhomo- 
geneous chemical potential. We argue that similar con- 
finement/deconfinment competition of intra- and inter- 
ladder interactions is at play. We would like to demon- 
strates this competition of interactions and the core effect 
where a spin gap and pairing is established on a sub-part 
of the system and it's spin degrees of freedom decouple 
from the rest of the system which may remain gapless. 
Such an effect is unique to models of coupled inequiva- 
lent ID systems. Below we examine the simplest model 
of such type. 

Consider a model of two inequivalent ID chains (or 
two inequivalent ID bands): Band- A has a Fermi wave 
number kA and attractive interactions —Ua- Band-B 
has a different Fermi wave number fc^ and repulsive 
interactions -hUs- The filling fractions are such that 
kA < kB < Hence, in the absence of inter-band in- 
teractions, the groundstate of band-A is a Luther-Emery 
(LE) liquid characterized by a spin gap (due to the at- 
tractive — Ua interactions) while band-B has gapless spin 
excitations. If band-B is half filled, fcs — ^, then it 
is equivalent to a Heisenberg spin chain, or otherwise if 
ks < ^ then band-B is a Luttinger liquid. Qualitatively, 
band-A and band-B correspond to the doped 2-legs and 
the remaining N — 2 legs of the inhomogeneous ladder 
model. ( The fact that in a 2-leg ladder the spin gap is due 
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FIG. 3: Fixed points and renormalization group flows. 



H {x) - {x) + Ubp"" [x) {x) 
+ (x) - UaP^ {x) / (x) 
+ Ms A [x] ■ SB (x) + J \p\ (x) Pb (x) + h.c. 



(1) 



pf (x) = V-A,! (^) ^A,^ i^)^ {x) = (a;) ^i^A,p (^)- 
Additional marginally irrelevant interactions (e.g., inter- 
band charge density interactions) which do not interest 
us here have been dropped. 

The three interaction {—Ua,M,J} are all marginally 
relevant in the initial stages of renormalization. Yet, the 
first crucial observation that we make is that the intra- 
band attractive interaction —Ua and the inter-band ex- 
change interaction M are competing, i.e., —Ua and M 
control distinct fixed points for which in the final stages 
of renormalization either M — > or Ua — *■ respec- 
tively. The fixed points and renormalization group flows 
are drawn in Fig-4. These results are derived using the 
bosonization method as summarized below. 

The fixed points of our model cam be inferred by obser- 
vation from the bosonized representation of the model ; 
Using the bosonized form of chiral (left and right moving) 
fermion fields ^1/^ {x) = F^y^ : e+'v^[«-(^)±<A-(^)l :, 



where da{x) = J^^dx'Il„{x'), and [na{x'), (paix)] = 
—iS{x' — x), a ^tii- The anticommuting Klein factors, 
{Fg., Frj'} = (5o-.cr', are needed for the proper anticommu- 
tation of fermions with different spin. (Normal ordering 
eliminates the short-distance regularization factor l/Va 
that is common in the literature). In the following, all 
operators are implicitly assumed to be normal ordered. 
Also, we absorb the infra-red cut-off p, factors into the in- 
teraction coefficients). As commonly done, we re-express 
the operators in terms of bosonic spin fields (psix) — 
— and charge fields (j>c{x) = -^[4'^ + ^ilj and 
correspondingly defined momenta and He- e.g., pair 
creation ^ Ig+iv^Sc (,Qg (^y^27r0s) and spin flip oper- 
ator s+ ^ lg+4V27res j,Qg (y/27r0s). The resulting form of 
the model Hamiltonian density is 



H (x) = (x) - ^ cos (\/8^0f ) 



to repulsive interactions while in this two-chain model the 
spin gap in band-A is due to attractive interactions is not 
of conceptual significance) . 

It is straight forward to determine that the only 
marginally relevant inter-band interactions (with respect 
to the non-interacting fixed point) are inter-band spin 
exchange Msa (x) ■ sb (x) and Josephson coupling (i.e., 

inter-band pair hopping) J P\ (x) Pb (x) + h.c. (where 

the j-band pairing operator Pj{x) = i'j^i^] i)- Note, in where 
the case that band-B is half filled then the inter-band 
Josephson coupling is perturbatively irrelevant. We shall 
discuss both kB < and fcs = cases. 

The Hamiltonian density which describe our one di- 
mensional (ID) two-band model is 



(2) 



(x) + ^ cos (y8^</.f ) 



M 

J 



cos 



h.c. + {marginally irrelevant interactions} 



Hi {x) = - 
1 

+ 2 



Ki {Hit 



(3) 



The difference between attractive and repulsive intra- 
band interaction is expressed by the Luttinger parame- 
ters > 1 and > 1 for the charge sector. It is use- 
ful to define 0± = ^ (0f ± 0^) and df = ^ {df ± 9^) . 
For the purpose of delineating the potential spin gap 
fixed points, one needs to keep track only of interac- 
tion terms which may develop a non-zero expectation 
value. Since by general theorems it always holds that 
(cos (\/47r07) cos (-\/47r(/)j)) = 0, then we need to keep 
only the cos [^/inO^) cos U/ATTcjjf) part of the exchange 
interaction. Also, since > 1 then the Ub interaction 
is irrelevant and can also be dropped. Hence we arrive 
at 



H = H^ + [x) - % cos (V4^ 



(4) 



2M 



■ cos 



(^^e-^ cos (%/4^0+) 



-I- ^ cos (^V47r6'c ) cos (^\/47r0+^ + cos (^\/47r(/)^ 

4- {irrelevant and marginally irrelevant interaction} 

The competition between M and Ua is inferred from 
the fact that at a fixed point where the inter-band in- 
teraction M is relevant then (cos (\/47r0j)) ^ and 
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(cos (-s/iTTt/ij)) = and hence the intra-band Ua interac- 
tion is irrelevant. The resulting phase diagram depicted 
in Fig. 3 shows the fixed points and renormalization group 
(RG) flows. In an expanded following publication we 
shall elaborate on the complete RG equations and esti- 
mates of the spin gap magnitude renormalization. 

As discussed by jj| , a fixed point is specified by the gap- 
less order parameters. The enumeration of gapless modes 
is straight forward by observation. Only the (j)^ charge 
mode remains gapless, as expected from an over all in- 
commensurately filled system. To distinguish between 
fixed points with the same number of gapless modes, it 
suffice here to focus on the pairing order parameters 



in part of the ladder system while diminishing spin ex- 
change coupling with the rest of the ladder system) leads 
at intermediate temperatures to the establishment of a 
spin gap on a sub-part of the ladder (the two hole rich 
legs), where the temperature scale and the magnitude of 
the spin gap are characterized by those of a doped 2- 
leg ladder, and are independent of the number of legs of 
the full ladder, (c) Consequently, at intermediate tem- 
peratures the ladder separates into two subsystems of 
which the spin dynamics is respectively decoupled, (d) 
At low temperatures, global coherence is established by 
the proximity effect (Josephson coupling) as expressed 
by the flow towards fixed point IV in the phase diagram. 



+ Lb^Rbi) + iRA]'LAi + La\Rai)] 
iV^e^ (,+\ (,- 



[{Rb]Lbi + Lb^Rbi) - {Ra^Lai + La^Rai 



Ate - 2 [{Rb]LaI 



Ato - 2 [{RbtLai 



sin {v^(j)t) sin (%/7r0; 
-^BT-^^i) + (Ra^Lbi + La]Rbi) 

{\^(t>t) sin {Vt^O; 
Lb]Rai) - {Ra]Lbi + La^Rbi) 
cos (^/¥6'; 



The subscripts e/o indicate even/odd parity under ex- 
change of band indices. 

Table-1 delineate the gapless modes at each fixed point. 
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If band-B is half filled, kB_= ^, then it is equivalent 
to a Heisenberg spin chain |y| and the fixed points are 
modified since band-B does not support gapless charge 
modes and the Josephson coupling J is perturbatively 
irrelevant. Hence, modes Ate and Ato are eliminated. 
For an elaborate discussion of FP-II in this case see lal . 



DISCUSSION 

From a general theoretical perspective, our modified 
ladder model has the following special properties: (a) 
The spin gap of the doped ladder model is higher than 
of the half filled ladder, (b) The predominant flow to- 
wards fixed point III (of increased pairing interactions 



We here argue that the model captures the core physics 
of the stripe structure in doped antiferromagnets. While 
the mechanism for the creation of stripes is not generally 
agreed upon, the end result is manifestly a state of local 
charge inhomogeneity of elongated domains of alternat- 
ing hole rich and hole poor strips, i.e., alternating strips 
with different effective local chemical potential. There- 
fore, a repeated arrangement of our inhomogeneous lad- 
der model is a model representation of the stripe struc- 
ture (where the width of the ladder corresponds to the 
charge period of the stripe) . Such a coupled ladder model 
has the following properties: At intermediate tempera- 
tures a spin gap is established independently on each hole 
rich stripe, while the hole poor strips may remain gapless. 
Residual interaction may lead to coupling between the 
hole poor regions across the spin gapped hole rich stripe, 
but at intermediate temperatures such coupling may just 
as well be too weak and leave the hole poor regions decou- 
pled and thus manifest only local AFM correlations. At 
low temperature, global coherence is established by co- 
herent pair tunneling, leading to a superconducting state. 

These characteristics of our model realize the goals of 
the general framework first presented in 0. Yet, in terms 
of the pairing mechanism itself, the core effect presented 
here is orthogonal to that in . The pairing mechanism 
proposed in (3] had the pairing arising out of the coherent 
interaction between hole rich and hole poor regions. In 
contrast, the core mechanism in our new model is a de- 
coupling of the dynamics of the hole rich and hole poor 
regions. The proximity effect introduced by comes 
into operation in our model only at low temperatures to 
establish global coherence (generating the flow towards 
fixed point IV) . Elsewhere we shall expand on the signif- 
icant advantages (compared with previous stripe based 
proposals) of the model and mechanism introduced here 
to address phenomenological properties of pairing and 
superconductivity in the cuperates. 

We thank Steve Kivclson and Jiang-ping Hu for valu- 
able discussions, and the support of NSF grant #4- 
444025-KI-21420-2. 
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